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ABSTRACT 
A  sound  performance  evaluation  methodology  is  essential  to  the 
design  and  use  of  computer  systems.  Among  performance  evaluation 
methods,  queuing  network  models  have  been  widely  recognized  as  cost  and 
time  effective  tools,  capable  of  predicting  ability  and  necessary 
complements  of  simulations.  Algorithms  for  efficiently  estimating  the 
throughput,  response  time  and  other  measures  of  performance  have  been 
devised  for  a  broad  range  of  queuing  network  models.  Although  this 
range  is  very  rich  and  adequate  for  modelling  many  realistic  features 
of  a  computer  system,  nevertheless  it  has  major  limitations.  An 
important  one  is  the  inability  to  model  pre-emptive  priority  scheduling 
disciplines.  Although  exact  closed-form  solutions  for  such  network 
models  are  open  reseach  problems,  several  approximation  techniques  have 
been  proposed.  An  important  one  was  Sevcik's  [Sevc,  77]  "shadow  CPU" 
technique,  *ich  admits  a  product  form  solution  and  gives  performance 
bounds  which  are  quite  tight  in  :Tany  circumstances.  In  this  paper,  we 
extend  the  shadow  CPU  model  to  an  optimal  approximation  of  the  priority 
model,  within  the  restrictions  of  the  product  form.  We  propose  a 
computationally  efficient  sub-optimal  approximation  technique  which 
gives  tighter  performance  bounds  compared  to  the  original  shadow  CPU 
technique.  We     also      investigate      the      limitations      of      such     shadow 

server^product   form   techniques   for  modelling   priorities. 
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1.    Introduct:  ion 

With  Che  advancemenC  of  technology,  computer  systems  have  become 
exceedingly  complex  and  performance  evaluation  can  no  longer  be  done  on 
an  infomal  basis.  In  addition,  detailed  simulation  of  complex  systems 
is  very  expensive.  Fortunately,  analytic  models  based  on  queuing 
networks  have  been  shown  in  several  cases'  to  provide  predictive  ability 
comparable  to  that  of  extensive  simulations,  yet  more  efficient  and 
cost-effective  ( [Giam,  76],  [Rose,  /6I).  In  addition,  a  lot  of 
approximation  techniques  have  been  dievis'e'd  to  "model  situations  for 
which  cl':Bed-form  computationally  efficient  sblutions  are  not  known. 
(See,  for  example,  [Koba,  74],  [Geie,- ''-^5^,'  [^e'.'ibo"','  7  4]  ,  [Sa,Ch,  75], 
[Sh,3u,   77],    [Reis,   81],    etc.)  ' '^'- '    '•'       ■—^'■'^    •  '■ ' 

In  the  past  few  years,  the  early  work  on  analysis  of  queuing 
networks  (see  [Jack,  63],  [Go.Ne,  6/],  [Klei,  75],  [Klei,  76], 
[  Ba ,  Mu ,  73])  has  been  substantially  extended.  Buzen  ([Buze,  73]) 
presenteti  a  very  efficient  computational  method  for  deriving 
performance  raeasures  from  cLased  queuing  networks  with  exponential 
servers  and  a  single  customer  class  ("product  form,"  first  analyzed  by 
[Jack,   63].  Basket,        Chandy,        Muntz        and        Palacios        ([BOP,    75] 

demonstrated  that  certain  more  general  networks  have  the  pr^^duct  form 
property  of  state  probabilities.  Although  the  class  of  queuing  network 
models  admitting  a  product  form  solution  is  very  rich  ([Court,  72], 
[Re,Ko,  7  5]),  nevertheless  it  has  major  limitations  of  which  an 
important  one  is  the  inability  to  accommodate  the  pre-emptive  priority 
scheduling  discipline  where  the  jobs  from  the  higher  priority  class 
pre-empt  the  jobs  from  tl-ie  lower  class  and  within  each  class  the 
scheduling  is   FCFS    (First   Cksme    First     Serve).        i^euing      networks      with 
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pue-emptive  prioricy  scheduling  caa  be  analyzed  by  Che  nie  Chod  of  global 
balance  equations  (see  [Kiel,  75]).  This  method  results  in  a  system  of 
linear  equations  which  can  be  numeri^cally  solved  to  determine  steady 
state  probabilities  of  the  system  states..  ..  l\owever ,  the  cardinality  of 
the  system  state  space  grows  exponejiCially  with  the  number  of  service 
centers,  customers   and  classes.  .  ._.    .^    _ 

Although  priority  schedulir\g.  is  very  imp,p-rte^nt  and  basic  in  most 
nul  ti-prograraming  computer  sys.te,ms ,  .no  e,5tact  qlo.sed  form  solutions  of 
networks  with  priorities  are  yet_kn,own.  Jhus,  r.esearchers  attempted  to 
efficiently  model  priority  scheduli;^  thr:<jugh.  approximation  techniques 
(see  [Mitra,  72],  [Reise,  76],  [Sa.Ch,  75],  [Sevc,  77]).  Sevcik's  work 
[Sevc,  77],  provided  a  sequence  of  .  aporq:.:iraa tions  using  queuing 
networks  with  product  form  state  probabilities.  Some  of  those 
approximations  used  the  idea  of  introducing  a  second  central  processor 
(the   shadow  CPU)   for  Che   exclusive   use   of   Che    low  priority   class. 

In  this  paper  we  extend  the  power  of  Che  shadow  CPU  approximaCion 
to  provide  an  optimal  approximation  to  Che  priority  model,  which  is 
pessimistic  for  the  performance  of  the  high  priority  class  and  admits  a 
product  form  solution.  We  also  provide  a  suboptiraal  model  which  is 
computationally  efficient  and  outperforms  the  original  (modified) 
shadow  CPU  of  [Sevc,  77].  We  also  investigate  the  limitations  of  Che 
technique  of  using  a  separate  server  for  the  lower  prioricy  class  and 
staying  within    the    product   form  class   of    networks. 
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2.    The   Priori  ty  Scheduling  Model  and  Sevcik.  's   Work 

We  shall  use  a  model  proposed  by  [Sevc,  77]  chac  contains  very 
lictle  of  Che  complexity  of  an  actual  syscs;,.  but  still  reflects  the 
intrinsic  difficulty  of  analytically  appros-chinj  priority  scheduling. 
The  queuing  network  is  a  closed  cycli.c  uet^-^ork  of  two  servers:  a  CPU 
and  an  I/O  server.  There  are  just  fwo  customer  classes,  with  class  1 
having  pre-emptive  priority  over  class  2_  at  the  CPU,  while  all 
scheduling  is  FCFS  at  the  I/O.  That  Is  ,..^JTe  central  processor  will  be 
devoted  to  class  1  whenever  a  class  1  job  1-s-  available  for  service. 
Within  each  class  jobs  are  served  in  •FCFS^xjjcder ., . -The  average  service 
times  at  the  I/O  are  the  same  (X^)  i:  r  Vjfbo<:h:  <fclas.ses-  but  the  average  CPU 
service  times  (between  I/O  requests)  are- .different  for  the  two  classes 
(X.  for  class  i,  i  =  1,2).  Jobs  always  reCurr  to  the  CPU  after  being 
served  at  the  I/O.  'fe  assume,  for  simplicity,  that  the  service  time 
distribution  of  the  I/O  is  exponential  (of  rate  U3  =  ^3  )  and  the  CPU 
service  for  class  1  is  exponential  (rate  y,  =  ]C  ).  For  class  2  the 
same  holds,  when  the  CPU  service  time  distribution  is  conditional  to 
the   event  of   absence   of   class    1    jobs    from   the    CPU. 

This  model  (called  the  exact  model)  has  the  essential  features  of 
a  computer  system  with  an  interactive  and  a  batch  class  of  jobs,  where 
the  interactive  jobs  have  pre-emptive  priority  at  the  CPU,  The  I/O 
server  is  a  crude  model  of  several  channel/device  groups.  Let  N.  be 
the  total  number  of  class  i  customers  (i  =  1,2)  and  let  N^  be  constant 
(fixed  multi-programming  level,  heavy  load  assumption).  Without  loss 
of  generality  we  will  quantify  the  quality  of  an  approximation  to  the 
exact  model  by  tlie  degree  it  approximates  the  average  response  time  of 
Che  class  1   jo'os  at   the  CPU. 
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[Sevc,  77]  proposed  that  tn  order  to  eliminate  contention  at  the 
CPU  due  to  class  2  jobs  (but  without  ignoring  them),  a  virtual  CPU 
should  be  introduced  (shadow  CPU).  The-  Shadow  CPU  exclusively  serves 
the  class  2  jobs  (with  mean  rate  U2''^~^  )•  ''^^  call  this  the  shadow 
CPU  (SCM)  model.  In  the  S  CM  model  -cTa^^^'i  -jobs  receive  unrealistical  ly 
good  service  at  their  own  CPU",' tHii^  SCM  is  a  pessimistic  model  for 
class  1.  Sevcik  proposed  that-  In  order  to'-get  a  better  approximation, 
we    should   increase   the   mean -«erv ice  %'tiiie  hi   th-e   class   2   jobs    from  X^    to 

1 — ,      where     p  ,      is     the"    CVM     utf-ilization      of      class      1.      By        tlus 

modification  (modified  ■  s-hadow  ^  flioVfel-',"-  MSMO  ^we  slow  down  the  CPU 
performance  of  class  2,  getting  ^a^'l^'ttfeV  rao<Je'l  for  class  I.  t^/hile  q^  is 
not  known  a  priori,  an  ite'ratl-ve^-aff^>tbac'ii-  could  be  used  to  reach  a 
self-consistent  class  1  CPU  XitllizaiC-ikjn.  MSM  differs  from  the  exact 
model  only  in  that  the  times  between  departures  of  class  2  jobs  are 
exponentially  distributed  over  the  busy  period  of  the  sliadow  CPU.  A.s 
[  Sevc,  77]  remarks,  in  the  exact  model,  the  class  2  interdepartur-e 
times  have  higher  variance  (since  they  are  interrupted  by  class  1  busy 
periods  at  the  CPU)  and  (since  increasing  variance  decreases 
throughput)  class  2  progresses  faster  in  the  MSM  than  in  the  exact 
model  and  produces  larger  I/O  queues,  causing  the  MSM  to  be  a 
pessimistic   model    for  class    1. 
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3.    The  Generalisation  of   the  MSM  Model 

To  provide  a  conceptaal  basis  of  the  various  results,  we  shall 
generalize  the  SCM  and  MSM  by  introducing  a  load  dependent  MSM  model 
(LDSM,  load  dependent  shadow  model)  in  which  the  mean  service  rate  of 
the      shadow      server,   when   i   class   2   customers   aie    in    the    shadow  server, 

is  U'^Ci)  =6(i)*U9  for  1  <  i  <  N2  with-  U2  ""  ^2'  ^^'^  ^^^^  ^  °  ^°^ 
1  <  i<  ^2.  The  functions  S(i)  (yet  to  be  specified)  allow  us  to  get  a 
spectrum  of  models  which  includes  the  pessimistic  as  well  as  the 
optimistic  models,  within  the  product-form  restriction.  (See  also 
[Sh,    79].) 

1-      zsc'~..-r      .    : 
Table   1^ 

(1)  all       B(i)    =0  =>       optimistic         model 

(2)  all       B(i)    =1  =>       SCM 

(3)  all       B  (i)    =    1t3  ^        =>       MSM 

(^)      0   <   3(i)    <    1~P  ,    =>     product    form   models    bounded 
by  the    pessimistic   MSM   and    the   optimistic  model. 

(5)      All   i,     L-p  j^  <   B(i)   <    1    =>  A.11   pessimistic  models 

between  SCM  and  MSM. 

>btice  that  the  exact  results  will  be  best  approximated  by  the  category 
(4). 


Lemma      1.        The  mean   service   rate   of   the   class   2    jobs    in    the   exact 


model   is  172   =U2(^    ~Pi)»    where  P  1    =  class   1    CPUutilizati 


on. 
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Proof.  1^2    ~  I    p(S)u7(S),      S  =  sEate      of      the      exact     model      and 

S 
UoCS)    =   0  if   S   implies    n^   >   0     where  ,.  nj. ,^,-1.3     the      number      of      class      1 


nee 


customers      at      the     CPU     and     U2(S^    =  ^2    ^^      ^"'  itnplies   n^    =  0.      He 
Mo    =       1  U2*P(S)  where  K   =    {S;   n^   =  P}  ,    i«e.  , 

s  e  K 

iTo    =  u-)*  Prob  (nj^    =  0)    =  U2(^    -  P  i) 


Let    us   denote  by    server  1,    2  or  3     the     CPU,      shadow      CPU     or     I/O 
respectivelv .  Let     a.     be      the      number     of     customers      at      server     j 

(j  =  1,2,3)  at  state  S .  Let  P^(i)  be  the  marginal  probability  of  i 
customers  at  server  k.  The  mean  rate  of  the  shadow  CPU  of  any  LDSM 
model   is 

^2 

U9   =     ^        P2(i)S(i)U7    • 
i=l 


Definition.      An  LOSM    is   called   reasonable   if   its  ^2   satisfies      the 


restriction: 


(^1)  P9   =  u^Ci    -Pi) 


'2   -  t'  2 


or,     in  other    terms 
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I      P2(i)0(i)    =  Pi(0)    =   1   -  Pi 
i=l 


For  example,  MSM  is  a  reasonable  LDSM  whereas  SCM  is  not.  Each  choice 
of  B=  {g(i),  i=  1,2,...,N2}  in  Table  1,  category  (4),  corresponds  Co 
a  produce  form  LDSI  whose  class  2  customers  still  depart  from  the 
shadow  server  independently  of  the  population  count  of  the  "real"  CPU. 
In  fact,  each  LDSM  has  an  "illegal"  (nonzero)  mean  transition  rate  of 
class  2  customer  departures  out  of  sta£"es  S-for  which  n,  >  0.  Let  us 
denote  the  probability  of  a  state  S  by  p(S,B)  to  indicate  that  n  is  a 
function  of  the  set  5  of  load  dependent  factors  of  the  shadow.  The  set 
of  feasible  system  states,  E ,  is  a  partition  of  two  disjoint  sets  E^ 
and  Z      where 


E+  =    {S  i£  E    with   n^    >   0}  ,        Z^  =    {S  C  E    with   n^    =   0} 


Let  E^  =    {S    il  E^  with   n2    =  i}  ,  E^   =    {S  c^'  E  ^  wi  th   n2    =  i } 


tor 


i  =  0,1,2, ...,N2,  and  E^  =  E^LE^.  Then  each  choice  of  3  implies  an 
illegal   mean   rate 

"^2 
(EQL)  ir  =     I        3(i)(        I  p(S,S)) 

For  any  set  B  of  positive  6(i)'s,  ir  is  positive,  hence  the  solution  of 
the  exact  model  cannot  be  obtained  by  any  choice  of  0(i)'s  in  Table  1. 
However,  the  exact  solution  can  be  approximated  in  a.i  optimal  way  if  we 
select  B  to  mi  nimi^.e  ir,  subject  to  the  constraint  of  reasonable 
LDSM's.      "fetice   that   the   constraint    can   be   i^itten   as 
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(R2)  V        S(i)    [        I  p(S,B)]    =         I  p(S,8) 

i=i  s  :  z^  s  -"^o 

i.e.  , 


^2  -    ."■■       ■^2 

Ir  +     I        B(t)    (        I  p(S,B))    =     I  I         p(S,B) 

o   _    ri  ''"■-     ^-^     --Ci=0  c     ,_r  i 


or 


^2 


(EQ 


2)  ir  =     I        (^    -S(i))    •    (        I""   p(S,B)]    +         I         p(S,B) 


i=i  •s  £^c*^-  ^      -    s.  z;, 


■-. :  n'-lz  1  T  :  e.c    . 

Corollary   1.      An    optimal   (with     respect      to     ir)      reasonable     LOSM 

approxiniatioa     to     the     exact  priority  model  can  be  obtained  by   solving 

the   following  rainimizatin  problem: 

Choose  ^  =  {6(i),  i=  1,  .  .  .  ,N2}  with  6(i)  >  0  to  minimize  Lr 
(given  by  (EQl))  under  the  restriction  (R2),  or  to  tninimze  ir  (given 
by   (EQD),    under  the   constraint    (EQl)   =    (EQ2)  . 

[3CMP,  75]  proved  that  the  probability  of  an  aggregate  system 
state  S   =    (npn2,n2)   is  of    the    form 

p(S,B)    =  C  •    gi(ni)   •    g2    (n2,B)    •    g3(n3) 

where      g^(n^  )      depends      only  on  the   conditions   at   server   i.    In   any  LDSM 
all  servers   are   of   type    1   hence    (e.g.,   for   the   shadow), 
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g2(n2,B)    =  eo"   


no         nj 
k=l 


where  eo  is  the  relative  visit  count  at  the  shadow,  (It  could  be  taken 
equal  to  Ut).  We  conclude  that  the  form  of  the  function  to  (uinimize  is 
(within   a    constant) 


A^  Ao  %2~^ 

BTTT  ^  6(1)6(2)   '^'  '  '^  S(l)6(2)...g-(it2-t)-      ..x  :   e ;■.:.. 

where  Aj^,  Ao,  ...  ,Aj^   _j^   depend   onNj^,    N2,  'u  j^ ,  'U  2  "^^  ^  3"       '^'-    ^'^      easy      Co 
prove    that 


Remark:      There      is    a   positive      i    =       r   (N,  jM-)  ,y  j^,y  2,ij  3)    (which   does 
not   depend  on  a   particular   state  S)   such    that    ir  <      ,   ir' ,    where 


••  r'    -       ^       +  ^  +       +  1 

Siry       6(1)6(2)      •"     6(l)...8(NO 


(To  prove  it,  consider  tlie  form  of  (EQl)  and  replace  numerators  (which 
are  <  1)  by  ones.)  We  shall  provide  a  suboptimal  solution,  which  is 
computationally  easy,    by  minimizing    ir'    instead  of   ir,   subject    to  Rl. 
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Theorem   U      The  set   3    =    (3  (i)  ,    i  =    1,  .  .  .  ,N2}   with 

2(^2  +1-1)  , 

6(i)    =  f .    . •     (1    -p,) 

N2(N2  +   1)  P2(i)  ^ 

.   .     bEni:     3'  :     ' 

rainimizes    ir'   with   respect    to  Rl..,--    _..  .    v-.  ,-      . 


Proof:  See  Appendix  1.  Note  that  Th'eqrenf  1-  gives  the  3(i)  as 
functions  of  the  unknown  quantities  P2(i)  and  p^.  .^  iteration  process 
is  needed  to  get  self -consistent  ,-8  Cii's  ^  .(of  .•  the  same  form  as  in 
[Sevc,    77]. 

Several      numerical      examples  showed      that    the 

(reasonable)  LDSfl  with  B  as  in  Theorem  1  is  a  pessimistic  approximation 
to  the  exact  solution,  better  than  the  MSf-I  under  most  load  conditions, 
and      of      the      same      computational      effort.  To        find        t'ne        optimal 

approximation  one  can  use  Lagrange  multipliers  and  Mewton-type 
numerical     solution      of      the      obtained      equations.  This        is        still 

computationally  more  efficient  than  the  global  balance  method,  since 
the  number  of  equations  is  proportional  to  M^  and  only  a  few  iterations 
are    usual Iv  needed. 


4.    Limitations    of    the  Product   Form  Shadow  CPU  Type  Models 

By  [BCMP,  75]  the  aggregate  state  probability  distribatlon  remains 
the  same  even  if  the  shadow  CPU  is  considered  to  be  a  type  2  server 
(instead  of  type  1),  since  the  state  descriptions  are  equivalent.  One 
Chen  tTiay  hope  that,  by  choosing  a  suitable  service  time  distribution 
for  Che  shadow  CPU  (with  rational  Laplace  transform)  we  would  be  able 
to  get  the  exact  solution.  This  is  unfortunately  not  so,  since  (for 
type  2  servers)  any  disCribuCions  (with  rational  transforms)  give  the 
same  results  as  the  exponential  service  time  distribution.  This 
implies  that  the  best  we  can  achieve  by  a  type  2  shadow  is  our  optimal 
approximation  .  model  of  Section  3.  It  is  impossible  to  get  the  e;<act 
solution  since  no  exponential  service  time  distribution  zeros  ir. 
However,  one  could  approximate  the  exact  solution  by  the  following 
iterative  scheme    (which  uses   PS   instead  of    FCFS    for   the    shadow): 

Process   ST^E 
initialization:      Let    the    shadow  CPU  be   as    in   the   MSM   model. 
Let    i  =    1 . 

1.  Find  the  1^  ^^,2"^,  .  .  .  ,i  +  1^^  moments  of  the  idle  periods  of  the 
class  1  server  (server  1)  based  on  the  solution  of  the  previous  model. 
Relate  this  to  the  l^'-,2"^,  .  .  .  ,i  +  1^^  moments  of  the  service  time  of 
the  shadow  CPU  by  inverting  and  multiplying  by  the  corresponding  moment 
of   an  exponential  variable   of   mean  X,. 

2.  Add  one  more  stage  (exponential)  to  the  shadow  CPU,  in  such  a 
way  to  obtain  the  moments  estimated  above.  (By  using,  e.q..  Cox's 
stages    representation). 
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3.  Go  Co  step  1  after  setting  i  =  i  +  I,  unless  you  have  agreement 
with   t'ue  "exact"   model    up   to  a    desirably  high   moment. 

The  networks  obtained  by  STAGE  are  product  form  (of  a  type  2 
shadow).  Hence,  by  our  previus  remarks,  they  can  never  be  out  of  the 
class  of  solutions  obtained  by  our  LDSM  approximations.  This  implies 
that  not  only  the  service  time  distribution  of  class  2  jobs  at  the 
exact  model  does  not  have  a  rational  Laplace  transform,  but  also  that 
its   higher  moments    are    important    in   determining    the    distribution.      This 

provides  evidence    that      the      exact      solution      is      not      a      product      form 

;  ;  '  '•  ir:  ^    •'.     -'■     . 
solution,    as  our  next    theorem  proves: 

■  ■  j  :ii  .0--    3cr->-j      sr  :      : 

Theorem  2:  The  exact  priority  model  does  not  admit  a  product  form 
solution,  even  if  no  assumptions  are  made  about  the  service  time 
distributions  of  the  CPU  and  I/O,  besides  the  pre-emptive  priority  of 
class   1    over  class   2. 


Proof.  We  shall  prove  that  the  exact  model  does  not  satisfy 
"Exponential     Local     Balance"      (ELB).      The   rest    follows   by  a   theorem  of 

[Ch,Sa,    78]   that   product   form   implies  ELB.      We  state    the   ELB   property: 

A  queuing  discipline   is   ELB    iff 

R(k/s)R(j/g_i^)    =  R(j/s)R(k/g_j) 
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for  all  s,  s-j,  s-k  states  in  Che  sCaCe  space  aad  all  classes  k  and  j, 
where  R(k/s)  is  Che  CoCal  raCe  of  service  for  class  k.  customers,  given 
s,    and   s-k   is   Che    sCaCe   wi  ch   one    less   class   k   cusComer. 

Indeed,  coasider  Che  state  s  with  one  class  1  customer  at  the  CPU 
and  an  I/O  queue  with  a  class  1  customer  first  and  a  class  2  customer 
second.  This  state  is  always  reachable  from  any  initial  state  of  the 
priority  model    provided   both  N,,N,  >    2»      We    have     R(.2/s'')    =   0     for      the 

total      rate   of    class   2   customers.      But  R(l/s''"')    >  0  and  R(2/s*  -   1)    >   0. 

*  .  .   ■'    -.---■■  ■*-  •  • 

Hence    for   s   =   s      and   k  =   2,    j    =    I    the  EL8— coreii  tioh    is    violated,    so    the 

,    ,      .  ,  r-  .  t  r  _—■-_.  c  .'       .- 

exact    model    t.s   not    product   form. 


a  Ci^: 
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appe:toix  i_ 

Proof  of  Theorem  J_: 
The  constraint   becomes 


^2 


I      e(i)  (^L_]    =  1 


By    the   fact   that   arithmetic   averages    of   positive   numbers   are    greater    tlian 
(or   equal   to)   geometric,   we   get 

1 

(I^J1)      ^  >    [0(1).[b(1)6(2)] (g(l)...S(N2))J      ^2 

2 

i-fence    ir'    gets    its  minimum  when  (INI)    becomes   an  equality,    i.e.  ,    when 

N2     N2-I 
6(1)   8(2)    ..  3  (No)   maximizes. 


From   the   form   of    the   constraint,    this   happens    when 


3(i)  P2Ci)     _  1  .-  11    , 

tor  all   i, 


N2  +  1   -  i)        1    "  P  1  "^2 

I  (N2  +   1   -  i) 

i  =  l 


2(N9  +   I   -  i)        1   -  pi 

=>  e(i)  = : •  1 

N2(N2  +   1)  P2(t) 


Note:    Given  a|^xj^   +  a9X2  +   ...  +  a>^x»^   =   1    (and  a^,x^   >  0)    the  product 
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■M 


Xm      ,   maximizes 


when 


^i^i 


=  same   for  all   I's   = 


1 


I      ^i 
i  =  l 
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